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Our ignorance of the dark energy is generally described by a two-parameter equation of state. 
In these approaches a particular ad hoc functional form is assumed, and only two independent 
parameters are incorporated. We propose a model-independent, multi-parameter approach to fitting 
the dark energy, and show that next-generation surveys will constrain the equation of state in three 
or more independent redshift bins to better than 10%. Future knowledge of the dark energy will 
surpass two numbers (e.g., [tuo,Wi] or [wo,w a ]), and we propose a more flexible approach to the 
analysis of present and future data. 



Standard candles such as Type la supernovae (SNc), 
as well as standard rulers such as the cosmic microwave 
background (CMB) and the baryon acoustic oscillation 
(BAO) scale, are currently the preferred probes of the 
expansion history of the Universe [l], @, H, Q] • By deter- 
mining distances at cosmological scales, these probes 
have firmly established that the expansion of the uni- 
verse is accelerating 0, S H B- It is now believed 
that a mysterious dark energy component with an en- 
ergy density ~70% of the total energy density of the 
universe is responsible for this accelerated expansion. 
The underlying physics of dark energy remains obscure 
and understanding the acceleration has become one 
of the foremost challenges in fundamental physics. 

In an attempt to discriminate observationally be- 
tween differing models of dark energy, it is useful 
to parameterize the dark energy by its equation of 
state (EOS), encapsulating the ratio of pressure to 
density. When model fitting data, it is generally as- 
sumed that the dark energy EOS follows a certain pre- 
determined evolutionary history with redshift, w(z). 
Common parameterizations include a linear variation, 
w{z) = wq + w z z or an evolution that asymptotes 
to a constant w at high redshift, w(a) = Wo+w a (l — a), 



with a the scale factor |ll|, [12 1 



Fitting data to an assumed functional form leads to 
possible biases in the determination of properties of 
the dark energy and its evolution, especially if the true 
behavior of the dark energy EOS differs significantly 
from the assumed form [l3j. The issues related to 
model-dependent studies of the dark energy EOS are 
a greater problem for the high precision datasets ex- 
pected from next-generation cosmological experiments, 
including distance measurements from a Joint Dark 
Energy Mission (JDEM [H]). 

Instead of using a parameterized form for w(z), one 
can utilize a variant of principal component analysis 
14[ to establish the EOS without relying on a specific 
parameter description of the underlying redshift evo- 
lution. This was applied by Huterer & Cooray [l5| 
to a set of early supernova data. More recently Riess 
ct al. Q used the same approach to analyze a new 



set of z > 1 SNe from the Hubble Space Telescope, 
while an analysis involving a larger combined dataset 
[Til l has been presented in Sullivan et al. [13] ■ Here 
we use this model-independent approach to study the 
extent to which future data will constrain dark energy. 
We find that more than two independent parameters of 
the EOS can be determined with next-generation sur- 
veys. Our results argue against claims in the literature 
that next-generation surveys can only determine two 
parameters of the EOS as a function of redshift [H| 
and we motivate a model-free approach to study dark 
energy. 

To encapsulate the range of possible future dark en- 
ergy surveys, we consider six different data scenarios 
and generate mock data for each one assuming a flat 
ACDM cosmological model Our datasets are: 

• Case A: A catalog of 200 SNe uniformly dis- 
tributed in redshift out to z = 1.8; in addi- 
tion, two BAO distance estimates at z = 0.2 and 
z = 0.35, with 6% and 4.7% uncertainties, respec- 
tively. This case approximates the current state 
of the data in SNe @, H, [l^] and BAOs @,[igj]. 

• Case B: A catalog of 300 SNe uniformly dis- 
tributed out to z = 0.1, as expected from ground- 
based low redshift samples, and an additional 
2,000 SNc uniformly distributed in the range 
0.1 < z < 1.8, as expected from JDEM or similar 
future surveys [37| . In addition to the two BAO 
distances described in Case A, five additional 
BAO constraints at z = [0.6,0.8,1.0,1.2,3.0] 
with corresponding fiducial survey precisions of 
[4.3,3.2,2.3,2.0,1.2]% (V1N1 from \2$). 

• Case C: The same SN dataset as described in 
case B, the seven BAO estimates as described 
in case B, and, in addition, ten new BAO con- 
straints expected from a proposed JDEM mis- 
sion by NASA or ESA concentrating primarily 
on BAO measurements, such as ADEPT (Ad- 
vanced Dark Energy Physics Telescope). These 
BAO estimates have precision (in Dy e 3j) of 
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[0.36, 0.33, 0.34, 0.33, 0.31, 0.33, 0.32, 0.35, 0.3^ 
0.37]% from z = 1.05 to 1.95 in steps of 0.05 

• Case D: A dataset of 10,000 SNe uniformly dis- 
tributed out to z = 2. In addition, seven BAO 
constraints as in case Rbut assuming stronger 
accuracies (V5N5 of [201) f° r the nve higher- 
redshift BAO constraints: [1.9, 1.5, 1.0, 0.9, 0.6]% 
at z = [0.6,0.8,1.0,1.2,3.0]. 

• Case E: The 10,000 SN dataset along with the 
seven BAO constraints as described in case D and 
the additional 10 BAO constraints expected from 
space-based missions as decribed in case C. 

• Case F: The SN dataset as described in case A 
combined with the BAO estimates in case E. 

For each of the cases listed above we create mock 
catalogs of SN and BAO observations. For each in- 
dividual SN we simulate a random distance modulus 
m{z) — M. consistent with our fiducial cosmological 
model. In cases A to C, we bin the Hubble diagram 
at z > 0.1 into 50 rcdshift bins, while the 10,000 SNe 
sample in cases A to F is binned into 500 bins. The 
error in distance modulus for each SN bin is given by 
a rn = V (< J int/%/AbiZ) 2 + S™?, where a- mt = 0.1 mag is 
the intrinsic error for each SN , N^- m is the number 
of SNe in the redshift bin, and 5m is the irreducible 
systematic error. We take the systematic error to have 
the form 8m = 0.02(0. l/Az) 1 / 2 (1.7/z max )(l + z)/2.7, 
where z m ax is the rcdshift of the furthest SNe, and Az 
is the width of the relevant redshift bin. This is equiv- 
alent to the form in [23}. In addition, we include the ef- 
fects of gravitational lensing magnification 24 , 25|, 2(| , 
although the noise from lensing is expected to be small 
for large datasets due to sample averaging 27 , 28[ . We 
make use of the probability distribution function for 
lensing magnification from Wang, Holz, & Munshi [2^] 
at z > 0.6, while for SNe at z < 0.6 we approximate 
the lensing by a Gaussian distribution for magnification 
with dispersion 0.093z magnitudes, as given in Holz & 
Linder [28j. 

We analyze a set of 10 independent mock catalogs for 
the data cases B and C, to account for random varia- 
tions in the estimates of the EOS. When model fitting 
each of the mock samples, we follow [3] and marginal- 
ize over a prior in Q m h (0.213 ± 0.023 [3(|), a prior 
in H a (72 ± 8 km s" 1 Mpc" 1 J3l[), and a prior in 
the distance to the last scattering surface at z = 1089 
(R = 1.71 ± 0.03 [HI). In order to account for uncer- 
tainty in the absolute magnitude, we also marginalize 
over the nuisance parameter, M y with a uniform prior 
of -0.6 to 0.6 @. To explore the importance of our fidu- 
cial assumption of a flat universe, for cases C and E we 
also explore a curvature prior, with a ler uncertainty 
on n k of 0.0032 j^. 

We refer the reader to Sullivan ct al. [T3| for details 
of our approach. We make use of a modified version 



of the publicly available wzbinned code [381 ] , which an- 
alyzes observational cosmological data via a Markov 
Chain Monte Carlo (MCMC) likelihood approach to 
estimate w(z) in rcdshift bins. In this parameteriza- 
tion, the EOS is taken to be constant in each redshift 
bin, but can vary from bin to bin. We take a total of 
six bins between z = and z = 3 and fix the EOS 
at higher redshift to a constant value of — 1 out to the 
CMB at z = 1089. In addition, we also impose a prior 
on our furthest (6th) bin: —5 < wq < 0. This bin re- 
mains largely unconstrained by the data cases we have 
studied, and the prior facilitates convergence. In what 
follows we will generally omit this bin, as it is poorly 
constrained. The rcdshift intervals of the first five bins 
are listed in Table 1 and the intervals are chosen so 
that the error on w(zi) is spread evenly across all bins. 

The integration to higher redshifts causes correla- 
tions between lower- and higher-redshift Wi{z) bins, 
and this must be taken into account in the subsequent 
analysis. These correlations are encapsulated in the co- 
variance matrix, which can be generated by taking the 
average of the Markov chain, C = (ww T ) — (w) (w T ), 
where w captures estimates of Wi{z) as a vector. This 
covariance matrix is non-diagonal, with correlations be- 
tween adjacent bins that slowly decrease with increas- 
ing bin separation. 

Instead of restricting ourselves to correlated values 
of Wi(z), we follow Huterer & Cooray [IH and decorre- 
late the EOS estimates. This is achieved by changing 
the basis through an orthogonal matrix rotation that 
diagonalizes the inverse covariance matrix. The Fisher 
matrix F = C _1 is then F = O t AO where the matrix 
A is the diagonalized inverse covariance of the trans- 
formed bins. The uncorrelated parameters are then 
defined by the rotation performed by the orthogonal 
matrix: q = Ow. 

There is a freedom of choice in the orthogonal ma- 
trix used to perform this transformation. We follow the 
approach advocated in [l5T ] and write the weight trans- 
formation matrix as W = O t A? O, where the rows are 
summed such that the weights from each band add up 
to unity. This choice ensures we have mostly positive 
contributions across all bands, an intuitively pleasing 
result. We apply the transformation W to each link in 
the Markov chain to generate a set of independent, un- 
correlated measures of the probability distribution of 
the EOS in each bin as determined by the observables. 
We denote these uncorrelated bins as w = Ww. When 
discussing our results, we will generally refer to these 
uncorrelated estimates. 

The results of our analysis are summarized in Table[IJ 
For each of the observational scenarios we give the la- 
errors on the determination of w in 5 uncorrelated rcd- 
shift "bins". Since the redshift bins are decorrelated, 
they now leak into one another; the new window func- 
tions are shown in Fig. 1 where we illustrate the ex- 
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FIG. 1: Uncorrelated binned estimates of the EOS from a 
typical mock sample generated for Case C (top panel) and 
Case E (bottom panel). The error bars show la and 2a un- 
certainties with solid and dashed lines, respectively. At the 
bottom, in separate colors, we show the window function for 
each of the independent, decorrelated bins. The overlap in 
window functions represents the relative contribution from 
adjacent redshifts bins to each uncorrelated w(zi) estimate. 



pcctcd errors in w(zj) for two mock samples from Cases 
C and E. For Cases B and C we also list the relative 
dispersion of the binned w(z) errors determined by an- 
alyzing 10 independent datasets. For Case C we find 
that the first four bins have average errors of 0.067, 
0.056, 0.048, and 0.059. Analyzing several indepen- 
dent mock catalogs, we found the standard deviation 
to be [3.3,1.6,2.5,2.9]%, respectively for the first four 
bins, relative to the mean error. In the last column of 
Table [H we specify the number (Np) of independent 
EOS parameters that could be determined, according 
to our analysis, to an accuracy better than 10%. This 
is the same criteria used previously to argue that next- 
generation surveys can determine at most two param- 
eters of the EOS [3. 

Our results show that it is possible to determine four 
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TABLE I: 68% error in the value of w in the uncorrelated 
redshift bins assuming a flat universe. The z range lists 
the redshift ranges for the original bins, but decorrelating 
the covariance matrix results in a leakage across bins. We 
show this leakage in the window functions in Figure 1 for 
two mock samples from cases C and E. The last column 
is the number of independent EOS parameters that could 
be determined to an accuracy better than 10%, consistent 
with a prior study [3- For cases B and C, within brackets, 
we list the scatter relative to the mean error of Wi(z) using 
a moderate number of random datasets corresponding to 
the two cases. For cases C and E, we also show the errors 
for the case where we allow for variations in curvature with 
a prior on Qk expected from Planck. 



independent EOS parameters to an accuracy better 
than 10%, for the case of 2300 SNe coupled with high 
precision BAO measurements (e.g., Case C; Fig. 1). 
With cases D and E, two of these four parameters are 
determined with an accuracy at the level of 5%. Fi- 
nally, we note that the results of case F show that even 
with only 200 SNe combined with future BAO mea- 
surements, one can constrain three independent EOS 
parameters to around 10%. 

The results in Table 1 are for the case of a flat uni- 
verse, except for an additional analysis of cases C and 
E with the inclusion of a Planck prior on the curva- 
ture [33j. With varying curvature, we find that the 
errors on w(z) broaden by less than 3% relative to the 
errors with the flat universe assumption. Thus the as- 
sumption of flatness, although it very slightly improves 
our fits, does not alter the general conclusion that one 
can constrain more than three parameters of the dark 
energy EOS. 

Linder & Huterer [ijj have argued that future data 
(SNe, CMB, and weak lensing measurements) will lead 
to a determination of no more than two independent 
parameters of the EOS to better than 10%. They 
consider a principal components analysis of the EOS 
binned at redshift intervals of 0.05, and argue that 
only the first two components are determined to better 
than 10%. While the third and higher principal com- 
ponents are determined invidually with less accuracy, 
by combining multiple components additional indepen- 
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dent precise w{zi) estimates can be achieved. Our ap- 
proach utilizes much wider (and uneven) binning in 
redshift, thereby allowing for a more robust capture 
of dark energy properties, and thus naturally finding 
multiple independent parameters. In addition, while 
(l8j limited themselves to a Fisher matrix approach, 
we explicitly generate mock data sets, accounting for 
intrinsic scatter and systematic biases such as lensing. 
Even with the inclusion of additional systematic errors, 
we find that future data constrain more than two in- 
dependent parameters of the dark energy equation of 
state. 

Two-parameter dark energy fitting has been incor- 
porated into the figurc-of-merit (FoM) quantity advo- 
cated by the Dark Energy Task Force [3], which is 
based on a two-parameter model for the equation of 
state (w(z) = wq + W a (l — a)). This provides a con- 
venient criterion for the evaluation of next-generation 
dark energy surveys. As an alternative, more general 
FoMs, such as the ones discussed in Albrecht & Bern- 
stein [35| and Sullivan et al. [l7|, allow for more than 
two parameters. Since different FoMs highlight differ- 
ent aspects of the theory and the data, consideration 
of a range of FoMs is warranted. 

In summary, we find that next-generation dark en- 
ergy surveys will be able to constrain three or more 
independent parameters of the equation of state to an 
accuracy better than 10%. This is in contrast to recent 
claims in the literature, and convenventional wisdom, 
that two parameters are sufficient in dark energy anal- 
yses. As we enter an era of precision measurements, it 
is important to avoid prejudicing our results with ar- 
bitrary functional forms for the dark energy. We have 
thus proposed a model-independent, multi-parameter 
analysis procedure for fitting the dark energy equation 
of state, and have shown that precision measurements 
of the dark energy can be expected with the next gen- 
eration of surveys. 
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